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Abstract. Based on the ISO spectral catalogue of compact Hn regions by Peeters et al. (2001), we present a 
first analysis of the hydrogen recombination and atomic fine-structure lines originated in the ionized gas. The 
sample consists of 34 Hll regions located at galactocentric distances between RGai= and 15 kpc. The SWS 
Hi recombination lines between 2 and 8 fj,m are used to estimate the extinction law at these wavelengths for 14 
Hll regions. An extinction in the K band between and ~ 3 mag. has been derived. The fine-structure lines of 
N, O, Ne, S and Ar are detected in most of the sources. Most of these elements are observed in two different 
ionization stages probing a range in ionization potential up to 41 eV. The ISO data, by itself or combined 
with radio data taken from the literature, is used to derive the elemental abundances relative to hydrogen. The 
present data thus allow us to describe for each source its elemental abundance, its state of ionization and to 
constrain the properties of the ionizing star(s). The main results of this study are as follows. The ionization ratios 
Ar ++ /Ar+, N++/N+, S +3 /S ++ and Ne ++ /Ne+, which measure the degree of ionization and to first order, the 
hardness of the stellar radiation, seem to increase with RGai- These ionization ratios correlate well with each 
other, implying that the spectral hardening affects equally the full range of ionizing energies. A Galactocentric 
gradient of N/O (AlogN/O = -0.056 ± 0.009 dexkpc" 1 ) is observed in the sense of a decreasing abundance ratio 
with Real in agreement with previous studies. Abundance gradients for neon and argon are derived of the form 
AlogNe/H = -0.039 ± 0.007 dex kpc" 1 and AlogAr/H = -0.045 ± 0.011 dexkpc" 1 . These elemental gradients 
could be enlarged by the existing Galactic T c gradient. Adopting a T c gradient of approximately 330 Kkpc -1 , 
the slopes in the Ne/H and Ar/H gradients become —0.06 and —0.07 dex kpc -1 , respectively. Lower limits for 
the sulphur and oxygen abundances are derived. Nitrogen abundances are derived for 16 sources. 

Key words. ISM: abundances - ISM: dust, extinction - ISM: Hll regions - Galaxy: abundances - Infrared: ISM: 
lines - Atomic data 



1. Introduction 

The distribution of Galactic elemental abundances is cen- 
tral to study the chemical evolution of the Milky Way. 
The relative abundances of the elements are sensitive to 
the star formation history, the number of massive stars, 
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the relative yield of the elements, and the exchange of 
matter between the disk and the halo through infall or 
ejection. Measurements of elemental abundances through- 
out the Galactic disk thus provide vital inputs to model 
the formation and evolution of the Galaxy. 

Hll regions are prime targets to derive the present- 
day elemental abundances. These objects consist of gas 
which is ionized and heated by the radiation of massive 
stars. When the massive ionizing star is newly formed, it 
is still embedded in its natal molecular cloud and the re- 
sulting H II region is a young, bright and compact nebula. 
Measuring the elemental abundances of the ionized gas in 
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Hn regions allows us to probe the interstellar medium 
(ISM) in the vicinity of massive stars and to trace the 
present composition of the ISM in the Galaxy. 

Abundance determination methods based on infrared 
observations present clear advantages over optical stud- 
ies because the infrared fine-structure lines are not very 
sensitive to changes in the electron temperature and do 
not suffer from high extinction. Infrared measurements 
can thus trace the elemental abundances in the central 
regions of the Galaxy. The infrared is also the only wave- 
length regime to measure the N ++ ion, the dominant form 
of nitrogen in highly ionized H II regions. 

Previous determinations of the Galactic distribution 
of elemental abundances based on the study of Hii r 



gion s have been performed in the optical ( 5 haver et al 



198 3|; ffich fc Silkey, 199l| ; Peharveng et al., 2000|) an d in 
the infrared, based on IRAS ( Simpson fc Rubin, 199C) and 



Kuipe r Airb o rne Observatory (K A O) observations ( Lester 
et al., 1 1987|; [Rubin et al., 1988j; [Simpson et al ., 1995b| 



Afficrbach et al., 1997] ; [Rudolph et al., 199?] ). These stud- 
ies resulted in a series of firm results concerning the varia- 
tion of elemental abundances across the Galaxy: the elec- 
tron temperature (T c ) of the H n regions increases with 
Rcai; the abundance ratio N/O decreases with Rcai; the 
abundances of the heavy elements N, O, Ne, S, and Ar 
decrease with Real- Additional optical studies are based 
on photosphcric emission lines from B stars ( Fitzsimmons] 
et al., 1992|; Smartt & Rollcston, 1997; Gummcrsbach 
et al., : 998|; Rollcston et al., 2000|; |Smartt et al., 2001|) and 



on planetary nebulae (Maciel & Koppen, 1994; Maciel & 



Quirez i, 1999). Similar results and gradients have been de- 



rived towards other galaxies (e.g. Vila-Costas & Edmunds 
199 2j; flaritsky et al., 199^ ) 



The Infrared Space Observ atory (ISO) spectral cata- 
logue of compact Hii regions ( Peetcrs et al., 2001 , here- 
after Paper I) presents the combined Short Wavelength 
Spectrometer (SWS) and Long Wavelength Spectrometer 
(LWS) grating spectra from 2.3 to 196 /jm for 43 nebu- 
lae. The catalogue tabulates the fluxes for the hydrogen 
recombination and atomic fine-structure lines. A detailed 
explanation on the error of these line fluxes, which will 
be used in the present paper, is given in Sect. 5 of the 
catalogue. The spectral coverage gives access to nearly all 
the atomic fine-structure lines in the infrared range. Lines 
from the elements C, O, N, S, Ne, Si and Ar are present 
in most of the sources. This paper will concentrate on the 
lines emitted by the ionized gas in the Hii region. The 
[Cn], [Oi] and [Sin] lines, produced by ions with ioniza- 
tion potentials lower than 13.6 eV, are expected to be 
mostly emitted in the Photodissociation Region (PDR) 
and will be analyzed by Damour et al. (in prep). For some 
of the ions (0 ++ , Ne ++ and S ++ ) two lines are present 
providing, in principle, an estimate of the electron den- 
sity. All the elements except for oxygen (note that 0°, as 
mentioned above, is produced in the PDR) are observed 
in two ionization stages (see Fig.|]), which alleviates the 
problem of applying ionization correction factors for un- 
seen ions (especially for S, Ne and Ar). Finally, the range 



of ionization potential covered by this set of data (up to 
41 eV) allows us to examine the ionization state of the 
H II regions and to constrain the properties of the ionizing 
star(s). From the whole sample, 34 Hn regions present 
enough emission lines to determine their elemental abun- 
dances and ionization properties. This subsample covers 
the Galactic plane from the centre to a galactocentric dis- 
tance, Real, of 15 kpc. It is thus possible to investigate 
trends of the ionization conditions in Hn regions and of 
relative and absolute abundances across a large part of 
the Galactic disk. Atomic parameters are crucial in order 
to interpret correctly the strengths of the different lines. 
The latest transition probabilities and collisional strengths 
have been compiled (cf. Table |l]) and used throughout this 
paper. 

This paper is structured as follows. Sect. || describes 
the kinematic distances and radio properties of the H n re- 
gions; Sect. H discusses the problem of extinction; a discus- 
sion of the densities of the Hii regions is given in Sect. |[ 
Sect. | gives an outline of the methodology used to de- 
rive elemental abundances with a particular emphasis on 
the advantages/disadvantages of the present data; Sect. || 
presents the results for the ionization state and its varia- 
tion with galactocentric distance. The ionic and elemental 
abundances, and their variation with Real, are presented 
in Sect. [?]; Finally, Sect. ^ discusses and summarizes the 
results of this paper. 

Table 1. References for the transition probabilities A and 
the collisional stengths fl used in this paper. IP refers to 
the IRON Project. 



Species 



Reference 



Nil 

NIII 

OIII 

Nell 

Nelll 

S III 

SIV 

Aril 

ArUI 



Calavi'c ot al, (1007 )-4Jf> XXII) 
Lcnnon fc Burke j 1994 ) (IP II) 
Galavi'i et al, (1998a)_aP XXIX) 



Blum fc Pradlian 1199 2) 



Gal 



Lcnnon fc Burk e (1994) (IP II) 



Mcndoaa (4fl&3^ 



Saraph fc 'fiTTR (1994) (IP IV) 



et al, 



IP XXII) 



Galavi'e el al, (1997 )41^ XXII) 
Bullcr fc Zcippcn 11994 ) (IP V) 

figsal) 



fff Rrnmajjr 



laval &: Gupt; 



Mriulny -i ( 1'IS'U 



TTTTTTOF 



Saraph S SToix y (1999) (IP XXX) 



Mo,»lf>-/» I 1'IX'i) 



fcian &; Berrm 



gton (1995 ) (IP IX) 



M^rWa 1.. 7.p; rr ™ (1983) 



FTOOTI MM £) 



Ualavis et al. (lyyrfb) (IP XXXII) 



2. Distances and radio properties 

2.1. Kinematic distances 

The study of the elemental abundance variation from 
atomic fine-structure lines emitted by H n regions relies 
on proper determinations of the galactic distances to the 
Hn regions and in this sense, a detailed literature study 
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157.4 Jim 4.9(1) 
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121.7 (im 2.7(2) 



[01] 

63.1 Jim 2.1(4) 
145.5 (xm 3.2(3) 



[Sill] 
34.8 Jim 
1.8(3) 



[NIII] 
57.3 jim 1.4(3) 



[OHI] 

51.8 Jim 3.1(3) 
88.3 Jim 4,5(2) 



[Nell] 


[Nelll] 




12.8 Jim 


6.1(5) 


15.5 Jim 


1.8(5) 


36.0 Jim 


2.6(4) 



[SIII] [SIV] 
18.7 Jim 1.0(4) 

33.4^ 2.2(3; 10 - 5 ^ m 3 - 7 < 4 > 



[Aril] 
6.98 Jim 3.7(5) 



[Arlll] 
9.0jim 2.1(5) 



Fig. 1. The fine-structure atomic lines observed in the combined ISO SWS/LWS spectra of compact Hn regions are 
shown as a function of the ionization potential. The electron critical densities, indicated for every line in italic, are 
given in units of cm -3 and are expressed as a(b) = a x 10 b . The critical densi ties for the ions with ionizati on potential 
lower than 13.6 eV (emitted in the photodissociation region) are taken from Tielens fc Hollenbach (1985| ). The other 



critical densities were calculated using the latest atomic parameters (see Table |l|) . 



has been done in order to get the most accurate distances 
to the program sources. 

The determination of the distances to the H II regions 
of the ISO sample has been discussed in Paper I and the 
reader is referred to this paper for detailed information. 
Table || lists in columns 2 and 3 the adopted distances 
Real and R Q (including in the latter the near and far 
solar distances). 



optical depth for an ionized, optically thin gas at a tem- 
perature T c is given by: 

T h = T.(l - e- T ) ~ T c t . (2) 
The emission measure can be derived from r using: 

T _g 

" .„ N _i 35 . — [pCCm ],(3) 

3.235 x 10- V (§) 35 (sfc) 21 



EM 



2.2. R.idio properties 



where a v is a cor rection factor close to unity ( Mezger fe 



Henderson, 1967 ) 



Radio continuum measurements allow one to derive vari- 
ous properties of the H n regions such as the emission mea- 
sure (EM), rms electron density (n C!rms ) and the Lyman 
continuum photon flux (A^yc)- These parameters can be 
derived for spherical sources if the radio flux density and 
the source size are known. For non-spherical sources, the 
geometrical average of the major and minor axes can be 
taken as a representative size. 

In the Rayleigh- Jeans limit, the brightness tempera- 
ture, lb, can be estimated from: 

1n -23 2 c 

where S u is the integrated flux density in Jy, v is the fre- 
quency in Hz and Q is the source solid angle in sr. The 



To estimate the r ms electron density we used the for- 
malism explained by Wood fc Churchwell (1989 ): 



'.EA/[pccm- 6 ] 



As[pc 



(4) 



where As is the path length of the emitting region (taken 
as the source size). Finally, the photon flux of the Lyman 
continuum can be determined from: 



7V Lyc = 1.54 x 10 55 (n e , rms ) 



As 
pc 



(5) 



where a B =3.28x 10" 13 cm 3 s" 1 ( [Hummer fc Storey, 1987|) 
is the case B recombination coefficient of hydrogen to all 
levels n>2 (for T C =7500K). 
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Table 2. Kinematic distances and radio properties of the sample sources. 



Source 


Rgal 


R Q W 


A 




Size 


EM 


n e , 


m 

ms 


l0g7V L yc (i,) 


Ref. 




(kpc) 


(kpc) 


(cm) 


(mJy) 


(") 


(10 6 pc cm" 6 ) 


(10 3 


cm -3 ) 


(s- 1 ) 




IR 01045 


13.8 


7.0 


6 


290 


3 


14.5 


11.9 




47.9 


1,2 


IR 02219 


11.0 


3.3 


6 


25000 


35 


9.2 


4.0 




49.2 


3,4 


IR 10589 


9.5 


8.0 


















IR 11143 


9.7 


8.8 


















IR 12063 


9.3 


9.5 


















IR 12073 


10.1 


10.8 


















IR 12331 


6.9 


4.5 


















IR 15384 


6.4 


2.7(11.5) 


















IR 15502 


4.6 


6.0 (8.4) 


















IR 16128 


5.5 


3.7(11.4) 


















IR 17160 


3.0 


5.7(11.0) 


6 


1300 


12 


4.0 


3.5 


(2.5) 


48.4 (48.9) 


5 


IR 17221 


5.2 


3.4 (13.4) 


















IR 17279 


3.4 


5.1 (11.8) 


6 


161 


10 


0.7 


1.7 


(1.1) 


47.3 (48.1) 


5 


SgrC 


0.3 


8.2 (8.8) 


21 


6400 


120 


0.2 


0.2 


(0.2) 


49.3 (49.4) 


6 


IR 17455 


0.5 


8.0 (9.0) 


6 


1600 


14 


3.7 


2.6 


(2.4) 


48.7 (48.8) 


5,7 


IR 17591 


5.5 


3.0 (13.8) 


6 


1310 


8 


9.2 


8.9 


(4.1) 


47.8 (49.1) 


7 


IR 18032 


7.6 


1.0 (15.8) 


6 


640 


15 


1.3 


4.2 


(1.1) 


46.5 (48.9) 


7 


IR 18116 


4.3 


4.5 (12.0) 


6 


3860 


20 


4.3 


3.1 


(1.9) 


48.6 (49.5) 


7 


IR 18162 


6.6 


1.9 (14.7) 


6 


4 


5 


0.07 


1.3 


(0.5) 


44.9 (46.7) 


1 


IR 18317 


4.5 


4.9 (10.6) 


6 


1290 


13 


3.4 


3.3 


(2.3) 


48.2 (48.9) 


5 


IR 18434 


4.6 


5.7 (9.0) 


6 


2280 


7 


20.9 


10.4 


(8.3) 


48.6 (49.0) 


5 


IR 18469 


4.8 


5.3 (9.2) 


6 


102 


13 


0.3 


0.9 


(0.7) 


47.2 (47.7) 


5 


IR 18479 


7.5 


1.2 (13.1) 


6 


2030 


7 


18.6 


21.4 


(6.5) 


47.2 (49.3) 


5 


IR 18502 


4.7 


7.1 


6 


1090 


6 


13.6 


8.1 


(2.2) 


48.5 (50.8) 


7 


IR 19207 


6.1 


5.7 


21 


13700 


140 


0.3 


0.3 


(0.1) 


49.3 (51.8) 


8 


IR 19442 


7.6 


2.5 (5.8) 


6 


23 


1 


10.3 


29.2 


(19.2) 


45.9 (46.6) 


9 


IR 19598 


9.8 


8.5 


3.6 


3475 


6 


45.6 


13.6 




49.1 


10 


DR21 


8.6 


2.8 


6 


17500 


20 


19.6 


8.5 




48.9 


11 


IR 21190 


12.7 


8.9 


6 


906 


4 


25.4 


12.1 




48.6 


12 


IR 21270 


14.8 


11.3 


6 


563 


120 


0.02 


0.1 




48.6 


13 


IR 21306 


12.6 


8.3 


















IR 22308 


11.3 


5.5 


6 


554 


90 


0.03 


0.1 




48.0 


13 


IR 23030 


11.4 


5.2 


6 


1900 


40 


0.5 


0.7 




48.4 


14 


IR 23133 


11.7 


5.5 


6 


1048 


15 


2.1 


2.3 




48.2 


14 



w Value in brackets is the far solar distance. 

^ Two different values arc derived for sources within the solar circle because of the near-far solar distance ambiguity. Values in brackets are 

calculated using the far solar distance. 

REFERENCES: (1) McCutcheon ct al. (1991 ) (2) |Rudolph ct al. (199(j|) (3) Roclf; 
(1994) 
(1989) 



h Goss (1991) (4) Ticftrunk ct al. (1997) (5) 3ccker et al 



((I) Liszt fc Spikcr (1995 ) (7) Saray ct al. (1993|) (8~ Mchringer (1994 ) (9) Barsouy (1989[ ) (10) Afflcrbach ct al. (199f) (11) Roclfscma ct al 
saacman (1984) (13) T ich (1993) (14) Fich (1986) 



Flux densities and source sizes were taken from the lit- 
erature selecting observations made at frequencies larger 
than 5 GHz. This ensures that the emission is optically 
thin. In addition, care was taken to select observations 
with a spatial resolution of the order of a few arcseconds 
to avoid structures larger than 10" to be resolved out. 
Radio properties derived using the above equations are 
listed in Table || for those sources with published observa- 
tions satisfying the above requirements. 

3. Extinction 

Compact Hn regions are heavily embedded in dust and 
gas and suffer a dimming of their radiation by dust lo- 



cated either within the ionized gas or in the neutral fore- 
ground material. The term "extinction" usually refers to 
the dimming of starlight, for which scattering from the 
beam contributes as much as true absorption because the 
star is effectively a point source. For extended objects, 
the scattering is less effective in attenuating the radiation 
because photons can be scattered into the beam as well 
as out of it. The effectiveness of scattering depends upon 
the geometry of the sources relative to the dust. In our 
objects, this is unknown. Hence, we will refer to the dim- 
ming of our objects due to dust absorption as "extinction" 
throughout this paper. 

The extinction towards compact Hn regions is gener- 
ally high with values of Ay ranging from ~ 10 to hundreds 



N.L. Martin-Hernandez et al.: Ionization and elemental abundances 



5 



of magnitudes - corresponding to extinction of ~ 1 to a 
few 10 mag. at near- and mid-infrared wavelengths. The 
extinction in the far- infrared is generally negligible. For a 
numbe- of Hii regions, the ISO SWS spectra show a se- 
ries of Hi recombination lines which can be used to derive 
the extinction curve characterizing the dust absorption to- 
wards 1 hese sources. Radio continuum emission can also be 
used tc derive the expected flux at any H I recombination 
line to provide an alternative estimate of the extinction 
when comparing with the observed line flux. 



Table 3. 

combinati 



Extinction Ak as determined using the H I re- 
lines method. 



Source A 



IR02 
IR10 
IR12 
IR12 
IR15 



K 



Source A 



K 



Source A 



K 



119 
389 
)63 
373 
384 



1.5 
1.5 
0.8 
0.8 
1.3 



IR 15502 
IR 17221 
IR 17455 
IR 18317 
IR 18434 



2.7 
1.8 
1.8 
2.0 
1.6 



IR 19598 
IR 21190 
IR 23030 
IR 23133 



1.6 
0.0 
0.0 
0.0 



3.1. Hi recombination lines method 

The mid-infrared H I recombination lines from the 
Brackett (n— >4), Pfund (n— >5) and Humpreys (n— >6) se- 
ries between 2 and 8 /jih can be used to derive the ex- 
tinction law at these wavelengths by comparing the ratio 
of observed line strengths to the predictions of recombina- 
tion theory assuming a nebular electron density, n Cl and 
temperature, T . The ratio of the hydrogen emission co- 
efficients is not very sensitive to n c or T c . Therefore, we 
have adopted rt c =1000 cm~ 3 and T c =7500 K, which are 
typical values for H II regions. In the following, we will con- 
sider the standard case B treatment, ie. where the nebula 
is optically thick to the Lyman series photo ns but opti- 
cally thin to all other series photons (see, e.g., Osterbrock, 



198S). The observed line flux, -F(A), can then be expressed 



as: 



^(A) = F (A)1(T - 4 ^ , 



(6) 



where F (X) is the non-extincted line flux, given by case B 
theory, and A\ is the extinction at the wavelength A. The 
relative extinction to a reference line (e.g. Bra at 4.05 /im) 
is: 



A A - A Bra = -2.5 log 



F(X) 
'F(Bra) 



. Fo(\) \ 
S F (Bva)J 



(7) 



The theoretical ratios are taken from Hummer & Storey 



(19871 ) 



Fig. H shows the relative extinction A\ — Aj3ra derived 
for the H n regions which have at least 3 H I recombination 
lines (14 sources in total) within the SWS range. The line 
fluxes are taken from Tables 7 & 8 in Paper I. 

At infrared wavelengths longward of 1.25 /zm, the ex- 
tinction law, as based on JHK photometry, appears to be 
independent of the line of sight and consistent with a sim 
pie power law A\ = A K (X/2.2)- 1 - 7 flMathis, 1990|; |Martin 



fe Whittet, 1990), were Ak is the extinction in the K band 
(2.2 /zm). For all the 14 sources we could derive an Ak 
for which this power law satisfies both the relative extinc- 
tion at 2.62 /jm (Br/3) and 4.05 /im (Bra). The solid lines 
in Fig. ^ show this power law for the corresponding Ak- 
Note that the curves generally fit the relative extinction 
at the other wavelengths nicely. The derived extinction 
values (Ak) are listed in Table 3. 

From the computed Ak one can derive the ex tinction 
at other wavelengths using the results tabulated in Mathis 



Table 4. Estimates of the extinctions A\/Ak to apply 
to the near- and mid-infrared lines using the "standard" 
extinction law described in the text. 



Line 


A( pirn) 


Aa/Ak 


Line 


A( (im) 


Aa/Ak 


Br/3 


2.6 


0.741 


[Neil] 


12.8 


0.165 


Bra 


4.0 


0.354 


[Ne in] 


15.5 


0.156 


[Aril] 


7.0 


0.140 


[Sm] 


18.7 


0.217 


[Ar in] 


9.0 


0.429 


[Sin] 


33.5 


0.069 


[Siv] 


10.5 


0.429 


[Nem] 


36.0 


0.060 



(1990), namely A g . 7 /A K = 0.544 for the 9.7 /im sili- 
cate feature and Ai$/Ak = 0.217 for the 18 /im silicate, 
in combination with the astronomical silicate profile of 
|Drainc (1985| ) and a A -2 power law for A > 20 /zm. Using 
this "standard" extinction law, the extinction Aa/Ak to 
apply to the near- and mid-infrared line fluxes is shown 
in Table g. In contrast to the extinction between 2 and 
7 /im, the extinction due to the silicate absorption at 9.7 



and 18 /zm depends on the line-of-sight (Draine 
and still both the shape and strength of these features re- 
main controversial. Therefore, the "standard" extinctions 
for wavelengths longer than 7 /im given in Table g arc 
merely illustrative and may not be applicable to all the 
Hii regions. 

Because of the lack of information on the extinction 
for most of the sample sources (extinction could only be 
derived for 14 sources out of 34) and the difficulties ex- 
plained above in correcting the lines beyond 7 /im, the 
line fluxes will not be corrected for extinction in the anal- 
ysis presented in this paper and instead the effect of the 
extinction will be indicated. 

3.2. Radio continuum method 



The intrinsic flux of any H I recombination line can be 
determined if both T c and the emission measure, EM, in 
the nebula are known. Radio continuum observations can 
provide this information (see Sect. ||). The intrinsic flux 
of Bra can be determined using: 



Fbtoi — ^BtoiEM [ — ) — , 
n, J 47r 



(8) 
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Fig. 2. Near- and mid-infrared extinction curves as derived from the Hi recombination lines. Each panel is labeled by 
the name of the Hn region and the derived extinction Ak- The solid lines show the extinction law (see Sect. |3.l| ). 



N.L. Martin-Hernandez et al.: Ionization and elemental abundances 



7 



-16.8 



5 



. -17.4 




IRAS 02219+6125 (W3 A) It is a very brig ht H II region lo- 
cated in a large 4' by 3' (Dickel et al., 1980) complex of H n 
regions, infrared sources and molecular clouds. Hay ward 



ct al. (1989) studied the extinction to W3A by compar- 



ing the 5 GHz map by Harris & Wynn- Williams (1976) 
with their infrared (H and K) observations. They found 
that the extinction varies from Ak ~ 1.0 to 1.7, in agree- 
ment with our computed extinction (Ak=1-5). The 9.7 
/im extinction was studied by Hackwell et al. (1978) ) , who 
found A9.7 to vary from ~2.0 at the position of W3 A to 
~4.0 at the perimeter of the W3 complex. This extinction 
due to the 9.7 /im silicate absorption is much larger than 
the Ag.7=0.8 derived from our Ak w hen the "standard" 
extinction law described in Sect. B.l is used. 



-18.0 



-17.5 -17.0 

iog(F B J - 



-16.5 



-16.0 



IRAS 12073— 6233 Near-infrared slit observations by 



radio 



Fig. 3. Comparison between the Bra flux predicted by 
radio observations and the extinction-corrected Bra flux 
observed by SWS. The solid line represents the one-to-one 
relation. 

where eera is the emissivity of Bra and £1 is the solid 
angle of the source. Using equations [l], 0, || and ||, the 
ratio of protons to electrons n v /n e ~ 0.90, T e =7500K 



and £Bra=l-39xlO 26 erg s 1 cm 3 ( [Hummer fc Storey 



1987j)7|we obtain: 



Armand ct al. (1996 ) mapped the variation of several Hi 
recombination lines across IRAS 12073—6233 and found 
Ak to vary from 1.2 to 1.7. We derive an extinction 
of Ak=0.8 mag. The discrepancy in ~ 0.4 — 0.9 mag. 
can be caused by the fact that we have observed the 
integrated H 1 emission in a 14" x 20" beam, which may 
lead to a lower hydrogen density column because of the 
contribution of extended emission around the nebula and 
thus, to the lower extinction we obtain. 

IRAS 18434-0242 (G29.96-0.02) It is one of the best 
studied ultracompact H11 regions and classified as 



cometary (Wood & Churchwell, 1989). By comparing im- 
ages of the 2.17 /im Br7 recombination line and 2 cm 



radio continuum emission (Fey et al., 1995), Watson 



where 'the flux densities; 3„, hum Table g ran be used. 

Fig. [| compares these predicted fluxes for Bra with 
the fluxes observed by SWS corrected using the extinc- 



et al. (1997) derived a map of the extinction towards this 



tion derived in Sect. ?7l: The radio flux densities and 



H 11 region. They found that the apparent extinction to 
the nebula is not uniform, with Ak varying fro m 2.2 to 
2.6 mag. Th ey give a mean value Ak=2.14±0.25. Pratap 



Abiq are considered to have a nominal 20% uncertainty. 
In general, the radio Bra fluxes are larger than the 
extinction-corrected SWS fluxes, especially for IR 02219 



et al. (1999 ) estimated Ak=2.2±0.1 towards the ioniz- 



ing star based on its H-K colour. As discussed above in 
the case of IRAS 12073—6233, the contribution of the ex- 



and IRj23030. These two sources are much more extended 
in radio than the SWS beam (see Table ||). On ihe other 

hand, the predicted radio line flux for IR 18434 is slightly 



tended envelope around the compact H11 region (Kim & 



Koo, 2001) to the SWS aperture may cause ISO to observe 



a lower hydrogen column density than the ones observed 
by the high ~1" resolution observations of |Watson et al 



of up to a few arc-minutes yKim fc Koo, 200l| ) which is — i ower a k =1 6 mag 



smaller This source Ifes In areglon of extended emission and [p ratap et aL (i^. This could explain our 



observed by SWS and it is not included in the quoted ra- 
dio flux density. This is probably the reason why the radio 
prediction leads to a lower line flux. 

The application of predictions from radio continuum 
observations to our ISO data depends on the source ge- 
ometry at both radio and infrared wavelengths and on the 
knowledge of the radio emission corresponding to the in- 
frared emission observed by ISO. We will return to this 
issue in Sect. ||. 

3.3. Comments on individual sources 

In the following, we comment on the extinction determina- 
tion in individual H 11 regions from which previous studies 
are available. 



IRAS 19598+3324 (K3-50A) It is the only Hn region in 
the sample for which the SWS spectrum was also taken 
at speed 4, yielding a resolving power of ~ 1500 (see 
Paper I) . In total, 20 H 1 recombination lines were detected 
(Tables 7 & 8 in Paper I) allowing us to derive a better 
sampled extinction curve in K3-50A than for the other 
nebulae of the catalogue, especially in the region from 2 
to 5 ^m (Fig. |^) . Note that the points around 3 /jm seem 
to be slightly above the standard extinction curve, which 
could be due to the 3 /jm O-H stretching mode seen in 
the spectrum (see Paper I). The extinction to K 3-50A 
has been studied in detail by |Howard et al. (1996| ) using 
both high-resolution near-infrared and radio images. They 
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found variations in the extinction across the nebula. They 
reported an extinction at 4.05 /im (Bra) varying from 
to 3 mag., with an averaged value of 1.6 in the central 
2.7" region. This is the same value we get from the Hi 
recombination lines analysis. 

4. Density 

The electron density (n c ) can be determined from the ra- 
tio of two atomic fine-structure lines of the same ionic 
species when they are emitted from levels with nearly the 
same excitation energy (e.g. Rubin et al., 1994[ ). This ra- 



Table 5. Electron densities derived from the [O III] 88/52 
/im line ratio. 



tio is sensitive to gas densities approximately in between 
the critical densities of each line (see Fig. |l|). A good in- 
dicator of gas with densities lower than a few thousands 
particles/cm 3 is the ratio of the [Om] lines at 52 and 
88 /im. Higher density gas can be probed using line ratios 
with higher critical densities such as [S in] 33/19 /im and 
[Nem] 36/15 /im. 

4.1. [Olll] diagnostic 

The densities derived from the [Om] line ratio for the 
ISO sample of Hn regions are in the range from ~ 100 
to 3000 cm' 3 (see Table ||). No dependence of the density 
on galactocentric distance is foun d (see Fig. |i|), in agree- 
ment with previous studies (e.g. Simpson et al., 1995b; 
Afflcrbach et al., 1997| ). 

The derived [Om] densities, n e ([Oni]), are compared 



Source 


n c ([OIII]) 


Source 


n c ([OIIIl) 




(cm" 3 ) 




(cm" 


) 


IR 02219 


2834 


-680 


IR 18116 


753 


+224 
-152 


IR 10589 


645 


+ 194 
-132 


IR 18317 


1543 


+738 
-397 


IR 11143 


290 


-71 


IR 18434 


817 


+260 
-171 


IR 12063 


1335 


^284 


IR 18469 


146 


+65 
-46 


IR 12073 


962 


+295 
-195 


IR 18479 


836 


+268 
-175 


IR 12331 


623 


+ 181 
-125 


IR 18502 


1181 


+411 
-257 


IR 15384 


1224 


+390 
-251 


IR 19207 


485 


+152 
-104 


IR 15502 


452 


+ 178 
-114 


IR 19598 


484 


+ 154 
-105 


IR 16128 


891 


+271 
-180 


DR21 


320 


+ 139 
-90 


IR 17160 


476 


+144 
-100 


IR21190 


2882 


+ 1987 
-871 


IR 17221 


594 


+ 182 
-124 


IR 21270 


186 


+74 
-52 


IR 17279 


127 


+61 
-44 


IR 21306 


214 


+ 109 

-70 


SgrC 


107 


+66 
-45 


IR 22308 


768 


+366 
-207 


IR 17455 


548 


+ 181 
-121 


IR 23030 


806 


+246 
-164 


IR 17591 


462 


+ 143 
-99 


IR 23133 


< 478 




IR 18032 


512 


+274 
-154 









calculated from the radio con- 



to the rms densities, n Cjrms , 
tinuum observations (see Table |J) in Fig. |S| In general, the 
rms densities are larger than the [O in] densities, except 
for three sources (IR 19207, IR 21270 and IR 22308), which 
are more extended (their radio sizes are 140, 120 and 90", 
respectively) than the LWS aperture. 

If the radio and ISO sample the same gas and there 
is no collisional de-excitation, n ([Om]) reflects the mean 
of the local density and must always be of the same or- 
der or greater than n Cjrms , which is derived assuming that 
the emission is spread uniformly along the line of sight. 
The ratio of both densities is the filling factor, which 
has been found to be typically of the or der of ~ 0.1 
for many Hn regions ( Copctti et al., 200C ). The sources 
with n c rms < n c ([Oin]) likely contain high density clumps 
embedded in a hot, low density gas from stellar winds. 
However, Fig. || shows that n c ([Om]) is much too low 
for most of the sources. At the densities implied by the 
radio continuum observations, the [O in] line at 88 /im 
(with n cr ; t = 450 cm~ 3 ) is collisionally de-excited in the 
ultracompact core of the Hn region. Considering that the 
real densities of these cores are probably even higher than 
our derived rms densities, even the [O in] line at 52 /im 
(iicrit = 3.1 x 10 3 cm -3 ) could also well be collisionaly 
de-excited. Therefore, the density derived from the [O in] 
lines is likely representative of a dilute shell surrounding 
the ultracompact core. This low density gas at high ex- 
citation (note that ++ is produced by ionizing photons 
with hv > 35 eV) could be produced by leakage of ra- 



diation from the core. Recent radio observations (Kurtz 
et al., 1999| ; Kim & Koo, 2001) have shown that the pres 



ence of physically related envelopes around ultracompact 
Hn regions is the usual case rather than the exception. 

4.2. [SHI] and [Nelll] diagnostics 

The derivation of electron densities using the [S in] and 
[Nem] line ratios from the ISO data yielded inconsistent 
results. Fig. || compares the observed line ratios of these 
ions with their expected values for uniform electron den- 
sities between 10 2 and 10 6 cm~ 3 , plotted as a solid line. 
The assumed T c of 7500 K has almost no effect on the 
plot. Meaningful [S in] and [Nem] line ratios should lie 
below their expected values at low densities, but we see 
that almost all points for the [Nem] line ratio, and many 
for [S in], lie above these limits by up to a factor of 3. This 
problem has also been rec ognized for the [Ne in] line ratio 
in some planetary nebula (Rubin et al., 2000). However, in 
that case the observed ratio is only 10% larger than what 
theory predicts for the low density limit. Possible causes 
of this discrepancy between observations and theory are: 

Aperture The lines involved in both [S in] and [Nelll] ra- 
tios have been observed in different apertures (14" x 27" 
and 20" x 33", respectively for the short and long wave- 
length lines). The black, large arrows in the bottom right 
corner of Fig. || indicate the shift in both line ratios that 
might be caused by aperture effect if the emitting region is 
larger than the SWS band 4 (20" x 33") and has a uniform 
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' 1000 



100 r 



and IR 19207), with a [S ill] ratio ~ 3, into the theoreti- 
cal limits. Instead, these sources are undoubtedly off be- 
cause only part of these sources were included in the SWS 
beam. Specifically, recent 6 cm ATCA maps of IR 16128, 
IR 17160 and IR 17279 (Martin-Hernandez et al. 2002, in 
prep) show that SWS was pointing at the edge of these 
sources and more of the source was included in the larger 
aperture of the longer wavelength lines than in that of 
the shorter wavele ngth lines. Likewise , radio maps of the 
sources IR 17591 flGaray ct al., 1993|), IR 1 8469 flKurtz 



et al., 1994[) and IR 19207 ( fvlchrmgcr, 1994|) reveal com 



5 10 15 

Rod ( k Pc) 

Fig. 4. Electron densities of the compact H n regions de- 
rived from the [Oin] 88/52 fj,m line ratio plotted against 
RGai- The arrow indicates an upper limit. 



1000 r 
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1000 
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Fig. 5. Comparison between the electron densities of the 
compact Hn regions derived from the [Om] 88/52 /im 
line ratio, n c ([0 in]), and the radio continuum, n Cjms . The 
dotted lines join the n c ^ ms as derived using both the near 
and far solar distances. The solid line represents the one- 
to-one relationship. The fact that n c ([Om]) is generally 
lower than n eirms shows that the [O III] lines are collision- 
ally de-excited in the central regions of the objects, so that 
ISO LWS observes only a surrounding shell of low density 
gas. 



brightness. Aperture effects may have affected the sources 
which cluster around a [S in] ratio of unity. Actually, this 
correction would bring these sources within the expected 
[Sin] and [Nem] ranges at densities of about 10 3 - 10 4 
cm -3 , the rms density range. However, this aperture cor- 
rection is not large enough to bring the most-offending 
sources (IR 16128, IR 17160, IR 17279, IR 17591, IR 18469 



plcx structure - core halo or multiple components - part of 
which was not included in the short wavelength, smaller 
aperture. 

Calibration The lines at the longest wavelengths, [S in] 33 
and [Ne in] 36 /im, are located in the SWS band 4. Band 
4 detectors suffer strongly from memory effects, are very 
sensitive to cosmic particle hits and their responsivity 
curve is still unsatisfactory and based on pre-fiight ground 
based tests (see Paper I). Therefore, these lines are af- 
fected by large calibration errors (an official la 25% cal- 
ibration error is quoted in Paper I and used in Fig. ^J). 
However, this calibration uncertainty is not enough to ex- 
plain the discrepancy between the observed and the pre- 
dicted line ratios. 

Atomic parameters The most recent theoretical calcu- 
lations (McLaughlin & Bell, 200C) of the collisional 



strengths for the [Ncm] levels differ from earlier studies 



(Butler & Zcippen, 1994) by only a few percent. Hence, 
we consider it unlikely that the observed factor 2-3 dis- 
crepancy in the [Nem] line ratios of Hn regions results 
from errors in the theoretical collisional strengths. 

Extinction The lines involved in both [Sin] and [Nem] 
ratios suffer also from differential extinction. The color 
excesses between the [S in] lines and the [Ne in] lines are 
E ~ 0.15 A K and E ~ 0.10 A K (cf. Table |), respectively. 
For the nominal Ak~ 2 mag. derived from the Hi recom- 
bination lines (cf. Table 3), these extinction corrections 
(indicated in Fig. ^ by the grey arrows in the bottom right 
corner) are, however, quite small. 

We conclude that the present [S in] and [Ne in] line ra- 
tio observations are unreliable as density indicators largely 
due to aperture differences. The sources above mentioned 
(IR 16128, IR 17160, IR 17279, IR 17591, IR 18469 and 
IR 19207), together with IR 21270 and IR 21306, which, al- 
though not present in Fig. ^| are known to be misspointed 
(see Paper I), will be excluded from the analysis involving 
any SWS line. Moreover, the [S III] 33 and [Ne III] 36 ^m 
lines will not be used in the data analysis. 

The object IR 21190, identified in Fig. ||, is of special 
interest. It is a very compact Hn region (and thus point- 
like for the SWS; see its spectrum in Fig. 8, Paper I, which 
does not present jumps between the different apertures) 
and has an Ak=0 (see Table 3). Thus, it should not be af- 
fected by aperture changes or extinction. Indeed, the [S in] 
and [Nem] ratios for this source are well inside the the- 
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.0 -0.5 0.0 
log([Nelll] 36.0/15.5) 

Fig. 6. Comparison between the observed line ratios [S in] 
33/19 /mi and [Nelll] 36/15 /im. Upper limits for the line 
ratios are indicated. The solid line represents the theoret- 



ical values of Lhe [3 III] and [NcIII] line i alius fui Lhe elee- 
tron densiLies 10 ? , 10 3 , 10\ 10 ri and 10 fi Tm . The laiu,e, 

black arrows in the bottom right corner indicate the shift 
(~ 0.23 in logarithmic scale) due to beam correction for 
an extended, uniform source. They are followed by short, 
grey arrows which indicate the extinction correction due 
to an Ak=2 mag. The source IR 21190 is indicated (see 
text). 



oretical limits. However, the densities derived from these 
two indicators (n e < 10 4 cm~ 3 and > 10 5 cm' 3 , respec- 
tively for the [S in] and [Ne in] ratios) do not agree. IRAS 
21190 presents a complex density morphology (Peeters et 
al. 2002, in prep) and it is therefore not surprising that the 
tracer of lower density and lower ionization state material 
yields a lower density than the tracer of higher density 
and higher ionization state material. 



5. Abundance determinations 

5.1. Methodology 

Ionic abundances can be determined from the measured 
strengths of the lines when a specific model for the struc- 
ture of the H II region is assumed. The simplest model is to 
assume that the nebula is homogeneous with constant T e 
and n c . If one also assumes that all the line photons emit- 
ted in the nebula escape without absorption and therefore 
without causing further upward transitions, the calcula- 
tion of the ionic abundances is straightforward (e.g. RubiiJ 
et al., 198$ ). 



where Fx+t and Fy+j are the fluxes corrected from extinc- 
tion corresponding to any line produced by the ions X +l 
and Y +J , and ex+ l an d ty+j are their respective emission 
coefficients. This equation assumes that the volume occu- 
pied by both ionic gases and the solid angles included in 
the fluxes are the same. Because the beam sizes of ISO dif- 
fer significantly with wavelength (see Paper I), the distri- 
bution of the emission over the beam becomes important 
in determining abundances for some of the elements. 

The emission coefficients depend on T c , n c and the rel- 
evant atomic parameters (transition probabilities, A, and 
collisional strengths, fi). They are computed using a 5- 
level atom (or a 2-level atom for ions with a 2 P ground 
term) in statistical equilibrium (Kafatos & Lynch, 1980). 
The references for the atomic parameters used in this pa- 
per are given in Table |l|. Fortunately, the emissivities of 
the fine-structure lines show only a very slight dependence 
on T c because these lines are emitted from levels with 
excitation energies much lower than the mean colliding 
electron energy. Thus, we will use an electron te mpera- 
ture of 7500 K, typical for galactic Hii re gions ( Shavei 



et al., 1983; Afiferbach et al., 1996, 1997), to evaluate 



these emissivities. However, the dependence of the fine- 
structure emission coefficients on n e is particularly impor- 
tant when the density of the emitting region is comparable 
to the critical density of the levels involved. It is illustra- 
tive to consider the n e correction factor, 8(X +t+1 / X +l ) nis , 
needed to apply to the abundance ratio of sequential ionic 
ionization states (Ne ++ /Ne + , calculated from [Neil] 12.8 
and [Nem] 15.5 fj,m, Ar++/Ar+, from [Aril] 7.0 and 
[Arm] 9.0 fim, S +3 /S ++ from [Sin] 18.7 and [Siv] 10.5 
fim, and N++/N+, from [N n] 122 and [N ill] 57 /im) when 
the respective line emissivities are evaluated in the low 
density limit (cf. Fig. ^). This correction factor is, ba- 
sically, the quotient between the ratio of the respective 
emissivities evaluated at a density n e and in the low den- 
sity limit. Over most of our parameter space (ie. n c < 10 4 
cm~ 3 , cf. Fig. ||), Nc ++ /Nc+ and Ar ++ /Ar+ can be con- 
sidered in the low density limit. The correction factor for 
S +3 /S ++ is not larger than 1.35. N ++ /N + can be overes- 
timated by up to a factor of 3. However, the strong depen- 
dence of N ++ /N + on n e can be circumvented by the use of 
the LWS [O in] densities, which characterize the ionic gas 
in the same region where the nitrogen lines are produced. 

To derive the ionic abundances with respect to hy- 
drogen, one needs the H + emission associated with the 
nebula. These can be derived from a near-infrared H I re- 
combination line, for instance Bra, which is detected in 
most of the sources, using equation (nw as follows: 



Fx+i/F Bla 

e X+ i / e BT a 



(11) 



Considering two ions X +l and Y" +J , the ratio of their 
ionic abundances is given by: 



X +i 



(10) 



The H I emission coefficients show a dependence on 
electron temperature tha t goes as T C ~ Q (a — 1.23 for Bra, 
|Hummcr fc Storey, 1987 ). The electron temperature in H II 
regions h as been shown to in c rease with the galactocentric 
distance ( {Shaver et al., 1983| ; |Afiferbach et al., 1996| , |1997 ; 
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log(n e [cm" 3 ]) 

Fig. 7. The n c correction factor S(X +i+1 / X +i ) n<! needed 
to apply to the ionic ratio X +t+1 / X +l when the line emis- 
sion coefficients are evaluated in the low density limit is 
plotted as a function of n e . 




log(n e [cm" ]) 



Fig. 8. The n c correction factor S(X +1 /H + )„ c needed to 
apply to the ionic abundance A +i /H + when the line emis- 
sion coefficients are evaluated in the low density limit is 
plotted as a function of n e . 



Dcharveng et al., 2000) from approximately 5000 K in the 
Galactic Center to 10000 K at 15 kpc. Hence, if a single 
electron temperature of 7500 K is assumed, the Galactic T e 
gradient will introduce a systematic bias in the abundance 
determinations. We can introduce a T e correction factor 
6(X+ i /R+) Te , defined as (A+7H+) To /(A+7H+) 7500K « 
(T c /7500)~ 123 . This correction factor goes from ~ 1.6 at 
5000 K to - 0.7 at 10000 K. 



The H + emission can also be derived from optically 
thin, radio continuum observations. Combining equations 
| and [ll]: 



X 



IT 



= 1.03xl0 _14 T " 35 



Sly 



m erg s 1 cm 3 , 



where Fx+* is in units of ergs cm , ex+i 
T e in K, v in GHz and S v in Jy. The ratio of electrons to 
protons, n e /n p , differs from unity because of the contribu- 
tion of singly ionized He to the total number of electrons. 
This can be included by using n c /np «l+N Ho+ /N H + w 
1.10. As can be seen from equation 12, this method shows 
a small dependence on T c and the assumption of a single 
T o =7500 K gives correction factors of the order of 15% at 
most. 

The dependence on n e of the fine-structure emission 
coefficients can also produce systematic uncertainties in 
the derived ionic abundances with respect to H + . Fig. || 
illustrates the case of the neon, argon and sulphur ionic 
abundances (calculated using the lines [Neil] 12.8, [Nem] 
15.5, [Aril] 7.0, [Arm] 9.0, [Sm] 18.7 and [Siv] 10.5 fj,m) 
and plots the n c correction factor, 8(X +% /H + )„ o sw 
[ex+'{n c )/ex+'(n c — ► 0)] _1 , needed to apply when the 
low density limit is assumed. This correction becomes 
larger than 1 when n e > n cr it- For the neon and argon 
abundances, this correction only becomes important for 
n c > 10 4 cm~ 3 . More critical is the case of the sulphur 
abundances because of the lower critical densities of the 
[S in] 18.7 ^m (10 4 cm" 3 ) and [Siv] 10.5 fim (3.7 x 10 4 
cm -3 ) lines. At densities around 10 4 cm -3 , the sulphur 
abundance will be underestimated by a factor of 2. Even 
more dramatic, although not shown in Fig. ||, is the case 
of N+/H+, N++/H+ and 0++/H+. As in the case of 
N ++ /N + , their strong n D dependence can be circumvented 
by the use of the LWS [O in] densities. 

The total abundance of a given element is calculated 
by adding the contributions of the different ionic species 
of that element. However, ionization correction factors 
(ICFs) must be applied when some of the ionic species 
are not observed. The nebular abundances are then deter- 
mined by: 



X 
H 



x- 



i>Q 



ICF uno bscrvcd ions X ^ 



X 



+> 



observed i 



(13) 



In the case of the ISO spectral data, all the elements 
except oxygen are observed in two ionization stages (see 
Fig. 1) and the ionization corrections will therefore be 
close to unity for Ar, Ne, S and N. In the case of oxygen, 
the abundance can only be derived by applying the ICF 
for the unseen + 

5.2. Photoionization models 

In order to evaluate in detail the dependence of the atomic 
fine-structure line flux ratios of the H n regions on the 
physical conditions (stellar and nebular properties), a 
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set of photoionization models based on the code NEBU 
( Morisset fe Pequignot, 1996| ) was performed. For the stel- 
lar spectral energy distribution (SED), the CoStar atmo- 
sphere models of schaerer & de Koter (1997) were used. 
These non-LTE stellar models include the effects of stel- 
lar wind and line blanketing, important for massive stars. 
A total of 26 CoStar models are available with effective 
temperatures ranging from 22 to 55 kK and luminosities 
from 1.5 x 10 38 to 8 x 10 39 erg s _1 . For each stellar atmo- 
sphere model, a series of photoionization models was build 
for nebular densities of 3 x 10 2 , 10 3 , 3 x 10 3 , 10 4 , 3 x 10 4 
and 10 5 cm -3 , which cover the density range probed by 
the atomic fine-structure lines detected in the ISO spectra. 
The nebulae are assumed to be spherically symmetric. Th e 
nebular abundances are solar ( Grevesse fe Sauval, 199g| ), 
the distance between the star and the nebula is taken to 
be 3 x 10 17 cm and the filling factor of the gas is supposed 
to be unity. In total, 156 nebular models were calculated 
spanning a wide range in stellar properties (one order of 
magnitude in luminosity) and nebular densities (two or- 
ders of magnitude), corresponding to a range in the ioniza- 
tion parameter of more than 3 orders of magnitude. This 
large range in ionization parameter combined with the ef- 
fective stellar temperatures is likely to describe the physi- 
cal conditions of the sample of H n regions. Since the main 
results deal with the predictions of line flux ratios, the 
distances to the objects are not relevant. A detailed com- 
parison of the predictions of such a model with the ISO 
spectrum of the compact Hn region IRAS 184 34-0242 
(G29.96-0.02) is given in pVIorisset et al. (2001| ) - here- 
after Paper III. 



6. Ionization state 

Any two successive stages of ionization i and i+1 of a 
given element can be used as a measure of the state of 
ionization of the nebula. The line ratios [Ann]/ [Aril] 
9.0/7.0, [N m] /[Nil] 57.3/121.7, [Siv]/[Siii] 10.5/18.7 and 
[Ne in] /[Neil] 15.5/12.8 /im can be used to determine 
(cf. equation |l^) the relative ionic abundance ratios 
Ar++/Ar+, N++/N+, S +3 /S++ and Nc++/Ne+, respec- 
tively. A single T e =7500K and the low density limit (ex- 
cept for N ++ /N + , which was evaluated using the [Olll] 
densities calculated in Sect. ||) are assumed in the deriva- 
tion. The results are presented in Table assuming no 
correction for beam differences or for differential extinc- 
tion. As mentioned in Sect. |^, the low density limit is 
valid for Ne ++ /Ne + and Ar ++ /Ar + over the range of 
the derived n e ([Olll]) and n o rms (n c < 10 4 cm -3 ), while 
S +3 /S ++ may be affected by up to a 35% uncertainty. The 
errors listed in Table [6] stem from the direct propagation of 
the line flux uncertainties (the N ++ /N + analysis includes 
also the uncertainties in the [Oin] density). [Ne ni]/[Ne n] 
is not affected by calibration errors because the lines in- 
volved in the ratio are measured in the same spectral band 
(see Paper I), resulting in smaller errors. 

Su ccessive stages of ionization X +l and X +l+1 of a 
given e lement indicate the state of ionization of the neb- 



Table 6. Ionic abundances X +l+1 /X + * for argon, nitro- 
gen, sulphur and neon derived in the low density limit 
(except for N++/N+, which was evaluated using the [O Hi] 
densities calculated in Sect. lil). 



Source 


Ar+ + /Ar+ 


N++/N+ 


S + 3 /s + + 


Nc ++ /Nc+ 


IR01045 


— 




< 0.07 


0.28 ± 0.07 


IR02219 


5.9 ±0.6 


> 0.8 


0.13 ±0.02 


1.07 ±0.06 


IR 10589 


0.90 ±0.10 


0.4±0.1 


0.018 ±0.004 


0.13 ± 0.01 


IR11143 


> 14 


3± 1 


0.21 ±0.04 


2.0±0.2 


IR 12063 


4.6 ±0.6 


> 1 


0.15 ±0.03 


0.88 ±0.06 


IR 12073 


15 ± 5 


> 3 


0.34 ±0.06 


2.7 ±0.3 


IR 12331 


2.1 ±0.3 


0.8 ±0.2 


0.027 ±0.005 


0.32 ±0.02 


IR 15384 


0.7±0.1 


0.5 ± 0.2 


0.020 ±0.004 


0.14 ±0.02 


IR 15502 


0.24 ±0.04 


0.5 ± 0.2 


0.009 ±0.002 


0.11 ±0.02 


IR 16128 


— 


0.7±0.2 


— 




IR 17160 


— 


0.31 ± 0.09 


— 




IR 17221 


0.28 ±0.04 


0.16 ±0.05 


0.004 ±0.001 


- ■ 


IR 17279 


— 


0.24 ±0.06 


— 




SgrC 


<0.04 


0.10 ±0.03 


<0.01 


<0.02 


IR 17455 


0.66 ±0.07 


0.6 ±0.2 


0.016 ± 0.003 


0.20 ±0.01 


IR 17591 




0.4±0.1 






IR 18032 


0.14 ± 0.02 


0.19 ± 0.08 


< 0.003 


< 0.01 


IR 18116 


0.12±0.01 


0.23 ±0.07 


< 0.002 


0.033 ±0.008 


IR 18317 


0.22 ±0.02 


0.10±0.04 


0.003 ±0.001 


0.029 ±0.002 


IR 18434 


0.56 ±0.08 


0.4±0.1 


0.018 ±0.004 


0.12 ±0.01 


IR 18469 




0.4±0.1 






IR 18479 


<0.1 


> 0.2 


<0.01 


0.21 ±0.03 


IR 18502 


0.38 ±0.06 


> 0.4 


0.011 ±0.003 


0.100 ±0.008 


IR 19207 




1.3±0.4 






IR 19442 


<0.05 




<0.03 


<0.01 


IR 19598 


1.4±0.2 


> 0.4 


0.09 ±0.02 


1.04 ±0.04 


DR21 


0.12 ±0.03 




< 0.009 


0.12 ±0.05 


IR21190 


2.9 ±0.4 




0.06 ±0.01 


0.43 ±0.02 


IR21270 




> 0.7 






IR 22308 


0.23 ±0.08 


> 0.2 


< 0.008 


<0.02 


IR 23030 


2.0±0.3 


0.5 ±0.2 


0.016 ±0.003 


0.063 ±0.003 


IR 23133 


0.25 ±0.03 




0.004 ±0.001 


<0.01 


ula, which 


depends 


only, in first approximation, on the 


ionization 


parameter, U , and the hardness of 


the ionizing 


radiation ( 


Vflchez & 


Pagel, 1988 


y 





oc 

v{X+ i 



^dv 



^dv 

13.6 cV hv uu 



(14) 



where J„ is the mean intensity of the radiation. Therefore, 
a ratio X +l+1 / X +l is, to first order and for a given U, pro- 
portional to the relative number of photons able to ionize 
X +t , ie. above v(X +r ) 1 as compared to that of Lyman 
continuum photons (above 13.6 eV). Hence, Ar ++ /Ar + , 
N++/N+, S +3 /S++ and Ne++/Ne+ probe the hardness of 
the stellar radiation between 13.6 eV and 28, 30, 35 and 
41 eV, respectively (see Fig. Q). 

The behaviour of these "hardness" indicators with 
galactocentric distance (cf. Fig. ^) shows that the high- 
est ionized H II regions are located at the highest galacto- 
centric distances. Assuming no systematic variations of U 
with Real, Fig. ^| indicates, consequently, a hardening of 
the radiation field towards large galactocentric distances. 

A ratio {X l+1 /X +i )l(Y 1+1 /Y +1 ) involving two differ- 
ent elements X and Y measures to first order, according 
to equation |lj, the relative number of ionizing photons 
in th e respective ionizing continua (e.g. Vflchez fe Pa; 



1988 ), ie. above v(X +l ) and v{Y +l ), respectively. Fig. 
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Fig. 9. The spectral "hardness" 
indicators Ar++/Ar+, N++/N+, 
S +3 /S++ and Ne++/Ne+ plotted 
against Real- The dotted lines are 
intended only to guide the eye: the 
vertical line indicates the solar po- 
sition at 8.5 kpc and the horizontal 
one indicates when X +l+1 / X +l = 
1. The upper and lower limits are 
indicated by down and upwards ar- 
rows, respectively. 



shows the correlations between the 4 "hardness" indica- 
tors, ie. log(Ar++/Ar+) vs. log(N++/N+), log(N++/N+) 
vs. log(S +3 /S++), log(S +3 /S++) vs. log(Ne++/Ne+) and 
log(Ne ++ /Ne + ) vs. log(Ar ++ /Ar + ). The arrows indicate 
the extinction correction to be applied if the "standard" 
extinction law (see Sect. |3.l[ ) is assumed. The extinction 
correction for the neon and nitrogen line ratios is negligi- 
ble. The dotted line in every panel represents a straight 
line fit to the data and is a fairly good representation of the 
data behaviour. The slopes obtained from a least squares 
fit are very close to one: 1.2±0.3, 0.76±0.07, 0.93±0.06 
and 0.8±0.1, respectively for panels a, b, c and d. Because 
these ionic pairs probe the ionizing energies up to 41 eV, 
Fig. |l^ shows that the hardening of the radiation field 
affects equally the full range of the ionizing spectra. 

Photoionization models, however, show an addi- 
tional dependence of (X +l+1 / X +l )/(Y +t+1 /Y +t ) on U 



( (5tasirisk"a~fc Schaerer, 1997 , Paper III) rendering the in- 
terpretation more complex and probably producing part 
of the scatter in Figs. ^ an d [To| . The set of photoionization 
models described in Sect. 5.2 allows us to study the influ- 



ence of the ionization parameter U and the stellar effective 
temperature T e g (ie. the shape of the ionizing spectrum) 
on the ionization state of the gas, traced by any of the ra- 
tios X +l+1 /X+\ for instance Ne ++ /Ne+. The variation of 
Ne ++ /Ne + is shown as a function of T e g in Fig. |ll|. For a 
given stellar atmosphere model or T e g, this ratio decreases 
with increasing electron density. The observed Ne ++ /Ne + 
ranges from ~ 0.01 to 3, a range predicted by the mod- 
els for values of T e g between 31 and 38 kK (Fig. |ll|) . 
Similar variations are predicted for the other ionization 



ratios (N++/N+, Ar++/Ar+ and S+ 3 /S++). Changing U 
at a given T e ff (by changing the stellar luminosity and/or 
the gas density) can change Ne ++ /Ne + by a factor of 10 or 
more. An estimate of the T e g of the ionizing star is there- 
fore only possible if robust constraints on U are available. 
A detailed study of the relations between the ionization di- 
agnostics and the stellar atmosphere models is postponed 
to a future paper (Morisset et al. in prep.). 

However, as tracers of the ionization state, the ionic 
abundance ratios X +t+1 / X +l are sensitive indicators of 
the ionization correction factors. Fig. [l2] shows the ICFs 
for N, O, Ne, S and Ar derived from the grid of photoion- 
ization models as a function of Ne ++ /Ne + . The scatter in 
the calculated ionization correction factors is low, indicat- 
ing that the use of the ionic abundance ratios X +l+1 / X +l 
to constrain the ICFs is a robust method. In other words, 
even if U and the SED of the ionizing star are poorly 
known, the knowledge of the ionization state of the gas 
through the ionization ratios is enough to determine reli- 
able ICFs. This conclusion does not strongly depend on 
the stellar atmosphere models as shown by comparing 
the results of a photoionization code using other atmo- 
sphere models (e.g. CMFGEN - Bourret, private commu- 
nication). 

The observed values of the ionization ratios can there- 
fore be used to correct the elemental abundances for the 
unseen ions. As shown in Fig. [l2], the effects are rather 
small, except for oxygen. The effect on the Ne abundance 
is negligible. In the case of sulphur, the abundance is al- 
ways underestimated by typically ~ 15%. For high ioniza- 
tion conditions, the abundance of N and Ar can be under- 



14 



N.L. Martin-Hernandez et al.: Ionization and elemental abundances 



: (a) 


































: -r 


1 





.0 -0.5 0.0 0.5 
log(N 4 "*/N" 1 ) 



; ( c ) 


* 






** 

















o -0.5 



: (b) 


/T " 












I •• t 


'4- 










► 



-2.5 -2.0 -1.5 -1.0 -0.5 0.0 0.5 
log(Ne 4 "*/Ne" 1 ) 




-1.0 -0.5 0.0 0.5 1.0 1.: 
log(Ar 4 "*/ Ar *) 



Fig. 10. Correlations between the 
ionic abundance ratios tracing the 
degree of ionization of the com- 
pact Hn regions: (a) Ar ++ /Ar + vs. 
N++/N+, (b) N++/N+ vs. S+7S++ 

(c) S +3 /S ++ vs. Nc++/Nc+ and 

(d) Ne++/Ne+ vs. Ar++/Ar+. 
Upper/lower limits are indicated 
by short arrows. The long arrows 
indicate the direction of extinction 
correction. If using the "standard" 
extinction quoted in Table |], the 
extinction corrections to apply are: (a) 
0.12xA K , (b) 0.08xA K , (c) 0.08xA K 
and (d) 0.12xA K . 



estimated by up to 10 and 30%, respectively, because of 
the missing N +3 and Ar +3 . For O, only ++ is available 
and thus, we are missing the low ionization state + . The 
ICF for O is very large, between 1.5 and 100. 

7. Abundances 

7.1. A/++/0++ and N/O 

The N ++ /0 ++ abundance ratio is straightforward to de- 
rive from equation 10 using the [Nm] 57 /Ltm and [Om] 
52+88 fj,m line fluxes together with the densities derived 
from the [O in] 88/52 line ratio (Table ||), which are well 
adapted because both the N and O lines are emitted by 
low density gas. All three lines are observed with LWS and 
consequently, beam differences and differential extinction 
are negligible. Finally, the use of both [O in] lines mini- 
mizes considerably the dependence of the ionic abundance 
on n e . The resulting N ++ /0 ++ ionic abundance ratios are 
listed in Table ^| (the quoted uncertainties come from the 
error propagation of the line fluxes) . The ionic abundance 
ratio N ++ /0 ++ is found to decrease with galactocentric 
distance (Fig. |l3| ). 

The significance of N ++ /0 ++ and its dependence on 
Real has been discussed in detail in previous studies (e.g. 
Rubin ct al., 1988, and references therein). N ++ /0 ++ is 



a function of both the ionization state of the nebula (the 
ionization potentials of N + and + are 29.6 and 35.1 eV, 
respectively) as well as of abundance, and one has to disen- 
tangle their effects to derive the N/O elemental abundance 
ratio. Fig. |lj presents the correlation between N ++ /0 
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Fig. 11. Model predictions for Ne ++ /Ne + as a function 
of stellar effective temperature (T D g) and nebular density. 
The densities are indicated by various symbols: the plus 
sign, star, diamond, triangle, square and circle correspond 
to densities of 3 x 10 2 , 10 3 , 3 x 10 3 , 10 4 , 3 x 10 4 and 10 5 
cm~ 3 , respectively. The observed range for Ne ++ /Ne + (~ 
0.01 to 3) in the ISO sample of Hil regions is delineated 
by the dotted horizontal lines. 

and Ne ++ /Ne + (note that the spectral "hardness" in- 
dicators are inter-correlated and thus, they all correlate 
with N++/0++). As can be seen in Figs. |l| and [ij, 
N ++ /0 ++ is inversely proportional to the degree of ion- 
ization (ie. low ionized nebulae have high N ++ /0 ++ ionic 
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Fig. 12. Model predictions for the ionization correction 
factors (ICFs) for N, O, Ne, S and Ar. The ICFs are plot- 
ted against Ne ++ /Nc + , which traces the ionization state. 
The dotted line indicates an ICF equal to one. The densi- 
ties are indicated by various symbols: 3 x 10 2 cm -3 (plus 
signs), 10 3 cm~ 3 (stars), 3 x 10 3 cm -3 (diamonds), 10 4 
cm~ 3 (triangles), 3 x 10 4 cm~ 3 (squares) and 10 5 cm~ 3 
(circles). 



Table 7. The ionic abundance ratio N ++ /0 ++ as derived 
from the LWS fine-structure lines [N in] 57 fim and [O in] 
52+88 fim. 



Source 


N++/0++ 


Source 


N++/0++ 


T"D aooi a 
IK Uzziy 


U.IO ± U.U4 


IK louoz 


n a _u n i 
U.O ± U.l 


TT3 1 nCiQO 

IK lUooy 


A 0/1 4- n C\A 

U.Z4 ± U.U4 


IK lollO 


n c^n 4- n as 
U.oU =h U.Uo 


T"D 1 1 1 A Q 


n i a _i_ a no 
U.14 ± U.Uz 


IK loot 1 


U.o ± U.l 


IK IzUoo 


a i o _i_ a no 
U.lo ± U.Uo 


TO A Q AO A 

IK 1o4o4 


a /i a _i_ a A'? 

u.4y ± U.U/ 


IR 12073 


0.13 ± 0.02 


IR 18469 


0.43 ± 0.07 


IR 12331 


0.30 ±0.04 


IR 18479 


0.15 ±0.04 


IR 15384 


0.41 ±0.06 


IR 18502 


0.39 ±0.07 


IR 15502 


0.47 ±0.08 


IR 19207 


0.24 ±0.04 


IR 16128 


0.38 ±0.06 


IR 19598 


0.10 ±0.02 


IR 17160 


0.50 ±0.07 


DR21 


<0.21 


IR 17221 


0.7 ±0.1 


IR 21190 


<0.33 


IR 17279 


0.7 ±0.1 


IR 21270 


0.17 ±0.03 


SgrC 


0.7 ±0.1 


IR 21306 


<0.13 


IR 17455 


0.44 ±0.07 


IR 22566 


0.2±0.1 


IR 17591 


0.12 ±0.02 


IR 23030 


0.33 ±0.05 


IR 18032 


0.6±0.1 







results, an empirical method can be determined to derive 
the elemental abundance ratio N/O from the observed 
[Nm] and [Om] line fluxes. The N/O abundance ratio 
can be expressed as follows: 



N 
O 



N++/0 H 



N/O 



[N in] 

S[Om] A [N iii]/E[Oiii] A N++/CF 



(15) 



where the first term on the right side is the ratio of the 
[N in] and [O in] line fluxes, the second term is the ratio 
of the emissivities of these lines (taking into account the 
density effect leading to collisional de-excitation) and the 
third term is the ratio of the ICFs which have to be applied 
to N and O. S[Oni] is the sum of the two [Om] lines, ie. 
[O m] 52+88 nm. 

The second and third member of the equation [l5| can 
be empirically determined by the use of appropriate ob- 
servables and the grid of models leading to the following 
relation: 



abundances) and this dependence is connected to the lo- n _ [Nm] 

O S[0 m 



cation in the Galaxy (ie. nebulae with high N ++ /0 ++ are 
preferentially located in the inner Galaxy). For these in- 
ner nebulae, N + and + are probably the dominant ionic 
species over N ++ and ++ and therefore, the ionization 
correction needed to get N/O from N ++ /0 ++ could be 
rather large (N ++ /0 ++ will be then overestimating N/O). 
In contrast, N ++ /0 ++ must be close to N/O for the neb- 
ulae in the outer Galaxy. 

The predictions of the photoionization models are use- 
ful to quantify the effects of ionization conditions on 
N ++ /0 ++ . In particular, one can estimate the ionization 
correction factor which has to be applied by using one of 
the ionization diagnostic X +l+1 / X +l . Based on the model 



x (l.07r[Oiii]-°- 29 
(0.95 + 0.271og 



JNcm] 

1 [Nc n] 



(16) 



where r[Om] is the ratio [Om] 88/52 /jm. The den- 
sity correction via r[0 ill] , which transforms [N m]/S[0 in] 
into N++/0++, yields N++/0++ ratios on average 10% 
lower than the N ++ /0 ++ derived from equation [l0[ 
mainly due to slight differences between the atomic con- 
stants used in the code NEBU and the ones used in this 
work. 

The photoionization model predictions of these suc- 
cessive ratios are shown in Fig. |l5| as a function of degree 
of ionization (as measured by the ionic abundance ratio 
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Fig. 13. Plot of the ionic abundance ratio N ++ /0 ++ ver- 
sus the galactocentric distance Real- The upper limits are 
denoted by downwards arrows. 

Ne++/Ne+): the line flux ratio [N m]/S[0 in] (plus signs), 
the ionic abundance ratio N ++ /0 ++ (diamonds) and the 
elemental abundance ratio N/O after correcting for the 
density and the unseen ions. The dotted horizontal line 
indicates the i nput value of N/O of the m odel, ie. the solar 
ratio of 0.123 ( Grevesse fc Sauval, 1998| ). As can be seen, 
when no ICF corrections are applied, the N ++ /0 ++ ratio 
overestimates the N/O abundance ratio for low ionization 
conditions and slightly underestimates the N/O ratio for 
high ionization conditions. Finally, the results for the N/O 
abundance ratio after correcting for the density and the 
ICF dependence (stars in Fig. |l5|) agree very well with the 
input solar value of 0.123. 

Equation [l6| provides a useful empirical law to de- 
rive the N/O abundance ratio from the observed [Oin], 
[Nm] and neon line fluxes. Any deviation from the solar 
N/O abundance after applying this law to the observed 
[Nlll]/S[0 in] indicates intrinsic variations in the N/O 
abundance ratio. Similar empirical laws can be derived 
for the [Arm] /[Aril] and [Siv]/[Sm] ionization diagnos- 
tics. However, we prefer the method based on the neon 
lines, whose ratio is hardly affected by differential extinc- 
tion and aperture differences. Note that the [N in] /[Nil] 
has an extra dependence on density and is, therefore, less 
appropriate. 

Fig. [l6] shows the N/O abundance ratio as a function 
of galactocentric distance after applying the corrections 
described above by using the neon lines. Note that there 
are only 25 sources for which the corrections could be ap- 
plied. Whereas the slope of the variation is not as steep as 
that for N++/0++ (see Fig. ||), a gradient of N/O with 
galactocentric distance is clearly present in the sense of a 
decreasing ratio with Rcai: 



N 

log- = (-0.23 ± 0.07) - (0.056 ± 0.009)R Ga i 
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Fig. 14. Plot of N++/0++ versus Ne ++ /Ne+, a degree of 
ionization indicator. 
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Fig. 15. Model predictions for the variation of 
[N ni]/S[0 in] (plus signs) as a function of the ion- 
ization diagnostic Ne ++ /Ne + oc [Nein]/[Nen]. The 
horizontal line indicates the input value of N/O of the 
model, ie. the solar ratio of 0.123. The ionic abundance 
ratios N ++ /0 ++ (after correcting the line fluxes for the 
collisional de-excitation) are shown as diamonds. The 
N/O elemental abundance ratios (after correcting for the 
density and the unseen ions) are shown as stars - see text. 



The correlation coefficient of the least square fit is 
r=— 0.80. This Galactic N/O gra dient is consistent with 
previous stu dies of Hii regions ( Simpson et al., 1995b ) 
and B stars ( Smartt et al., 2001 ). 



7.2. N/H andO ++ /H^ 



Because of the large difference between the SWS and LWS 
(17) apertures and the fact that many sources are more ex- 
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Fig. 16. Plot of the elemental abundance ratio N/O 
versus galactocentric distance Real as derived from the 
N ++ /0 ++ ratio - see text. The upper limits are denoted 
by downwards arrows. 

Table 8. Ionic and total abundances derived from the 
LWS fine-structure lines. These abundances are very un- 
certain because of the possible mismatch between the gas 
observed by LWS and the radio emission used to derive 



H+ (see Sect. 7.2) 



Source 


N + /H+ 


N+ + /H+ 


N/H 


O++/H+ 




xlO- 4 


xl0~ 4 


xl0~ 4 


xl0~ 4 


IR02219 


<0.3 


0.2 ±0.1 


<0.5 


1.4 ±0.6 


IR 17160 


1.5 ±0.6 


0.5 ±0.2 


2.0±0.8 


0.9 ±0.3 


IR 17279 


9. ±4. 


2.0±0.6 


11. ±4. 


3.0 ±0.9 


IR 17455 


1.3±0.6 


0.8±0.3 


2.1 ±0.8 


1.7±0.6 


IR 17591 


0.3±0.1 


0.12 ±0.04 


0.5±0.2 


1.0±0.3 


IR 18032 


1.7±0.9 


0.3±0.1 


2.±1. 


0.5 ±0.2 


IR18116 


0.5 ±0.2 


0.11 ± 0.04 


0.6±0.2 


0.21 ±0.07 


IR 18317 


4. ± 2. 


0.4±0.2 


5. ±2. 


0.5 ±0.2 


IR 18434 


1.3±0.6 


0.5 ±0.2 


1.9±0.8 


1.1 ±0.4 


IR 18469 


7. ±3. 


2.5±0.8 


9. ±4. 


6. ±2. 


IR 18479 


< 0.3 


0.07 ±0.03 


<0.4 


0.4±0.2 


IR 18502 


< 0.5 


0.20 ±0.08 


<0.7 


0.5 ±0.2 


IR 19598 


< 0.1 


0.05 ±0.02 


<0.2 


0.5 ±0.2 


DR21 


<0.04 


< 0.009 




0.04 ±0.01 


IR21190 


< 0.7 


< 0.2 




0.6±0.3 


IR 23030 


0.3±0.1 


0.16±0.06 


0.5±0.2 


0.5 ±0.2 



tended than the SWS aperture, the hydrogen line flux cor- 
responding to the gas where the nitrogen and oxygen lines 
are formed will be underestimated if a SWS H I recombina- 



cluded. [O III] densities (cf. Table |J) are used to evaluate 
the line emissivities involved. 



o 



O 




10 



15 



R 



Gal 



Fig. 17. (a) N/H and (b) Q++ /H+ as a function of R Ga i - 
The solar nitrogen abundance (|Grevesse fc Sauval, 1998 ) 
is indicated by the symbol 0. 



For nitrogen, both ions N + and N ++ are measured. 
Fig. [| shows that N + is the dominant species over N ++ 
and therefore, N +3 is not expected to contribute signif- 
icantly to the total abundance of nitrogen. Indeed, the 
grid of photoionization models (cf. Fig. [l2]) shows that 
the contribution of N +3 for the highest ionized sources 
is only at the 10% level. Hence, the total nitrogen abun- 
dance N/H is calculated by adding the ionic abundances 
N+/H+ and N++/H+. The calculated N+/H+, N++/H+ 
and N/H are given in Table |. Fig. |7| (a) represents N/H 
versus Real- An attempt to fit a straight line was made, 
but no statistically significant fit was possible because of 
the few points and the large scatter. We note that sources 
in the inner 5 kpc have supersolar nitrogen abundances, in 



tion lir e is used. Thus, the resulting abundances of N and 



O will be overestimated. Hence, the radio continuum mea- 
surements are thought to be more appropriate to evaluate 
H+ and equation [l2j is used to calculate N+/H+, N++/H+ 
and 0++/H+ via the lines [Nil] 122, [Nm] 57 and [Om] 
52+88 /xm. The radio flux densities from Table § are used 
and the sources which are more extended than the LWS 
aperture (SgrC, IR 19207, IR 21270 and IR 22308) are ex- 



agreement with determinations using B stars (e.g. 3martt 
et al., 2001| ) . 

For oxygen, only ++ is available and thus, we are 
missing + . We have used the grid of photoionization 
models to estimate the effect of the missing + abun- 
dance. These models show that the ionization correction 
factor is very large, between 1.5 and 100 depending on the 
degree of ionization (cf. Fig. |lj). Therefore, the derived 
++ /H + listed in Table || can only be interpreted as a 
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lower limit for the total O/H abundance. Fig. [l7] (b) plots 
++ /H + versus Real- Clearly, the derived ++ /H + abun- 
dances are much lower than the solar log(0/H)= —3.2 
( prevesse fc Sauval, 199g ). 

A note of caution must be brought concerning the ni- 
trogen and oxygen abundances. In most cases, the ISO 
[Oni] 88/52 /im line ratio give densities which are far be- 
low those required to produce the ultracompact H n emis- 
sion observed at radio wavelengths (n c , rms ) (see Fig. ||). 
This suggests that the ultracompact H II regions are sur- 
rounded by a tenuous shell of ++ (and probably also 
N ++ ). If this is the case, both N and ++ abundances 
will not be the true nebular abundances as the H + emis- 
sion associated with the shell will not be the one given by 
the radio observations. 



7.3. Ne/H, S/H and Ar/H 

Despite the inherent uncertainty associated with the H I 
recombination line method concerning T c (cf. Sect. |^), we 
will use the Bra line fluxes observed with ISO/SWS in an- 
alyzing all the SWS fine-structure lines of neon, sulphur 
and argon because the beams are well matched. As the 
lines of [Neil] 12.8 and [Nem] 15.5, [S in] 18.7 and [Siv] 
10.5, and [Aril] 7.0 and [Ann] 9.0, used to determine the 
abundances of neon, sulphur and argon, respectively, are 
observed by a beam comparable to Bra, the use of equa- 
ls straightforward. The fine-structure lines and H I 
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tion 

recombination line emissivities involved in the calculation 
were evaluated in the low density limit and considering 
a single T o =7500 K. The lines in the SWS band 4, [S III] 
33.4 and [Nem] 36.0 /jm, which are observed in a larger 
aperture and are affected by large calibration uncertain- 
ties, are not included in the determination (see discussion 
in Sect. |). 

The absolute abundances for neon, sulphur and argon 
were calculated adding the contribution from the different 
ionization stages available: 



Ne 

IT 
_s 

H 

Ar 

H~ 



Nc+ 

S++ 
H+ 



Ne^ 



H+ 
S+3 
H+ 5 



(18) 
(19) 



5 Lr + Ar 4 



H J 



H+ 



The grid of photoionization models described in 



Sect. 5.2 allows us to verify the validity of the above ex- 



pressions (cf. Fig. |12|) . The ionization correction factor for 
Ne is ~ 1 for all the sources. Ar +3 may be present in the 
higher ionized sources (photons with energies > 40.7 eV 
are necessary to ionize Ar ++ to Ar +3 , roughly the same 
energies to produce Nc ++ ). However, the contribution of 
this ion is estimated to be < 30%. The S abundance is 
always underestimated by ~ 15%. Table ^ gives the ionic 
and total elemental abundances for argon, sulphur and 
neon. The uncertainties are derived propagating the line 
flux errors. 



Fig. [l8] shows Ne/H and Ar/H as a function of R,Gai- 
Clear gradients with galactocentric distance can be ob- 
served for both species. Linear least squares fits to the 
neon and argon abundances, displayed as solid lines in 
Fig. are: 

Ne 

log— = (-3.49 ± 0.06) - (0.039 ± 0.007)R G ai and (21) 
H 

logTT = (-5-10 ± 0.09) - (0.045 ± 0.011)R Gal , (22) 
H 

with correlation coefficients r=— 0.76 and —0.66, respec- 
tively. 
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Fig. 18. Total abundances of (a) neon and (b) argon as a 
function of Real- The solar abundance ( jGrevesse fc Sauval 



1998) i s indicated by the sy mbol and the Orion abun- 



dance ( Esteban et al., 1998 ) by an open diamond 



We note that an electron temperature gradient in the 
Galaxy may systematically enlarge the above derived ele- 
mental abundance gradients. Specifically, if we correct for 
a T c gradient given by T c = 5000 + 5 000 x (R Gal (kpc)/15 ), 
based on the T c gradients de rived by |5haver et al. (1983 ), 
Afflcrbach ct al. (1996| ) and peharveng et al. (2000| ), the 
slopes in the Ne/H and Ar/H gradients become —0.06 and 
—0.07 dex kpc -1 , respectively. 

While Ar/H and Ne/H are Independent through the 
Bra line emissivity, the ratio Ar/Ne is practically insensi- 
tive to changes in T c . An approximately constant Ar/Ne 
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Fig. 19. S/H abundance as a function of Real- The solar 
abundance ( prevesse fe Sauval, 1998 ) is indicated by the 
sym bol and the interstellar abundance (Snow fc Wittj 



1996 ) Ipy a star. As discussed in Sect. 7.3, S/H could be 
underestimated by up to a factor 4 due to the T e and n e 
uncertainties. 



close to solar is found, contrary to the tentative Ar/Ne 
galactic gr adient derived from a s mall sample of KAO ob- 
servations ( Simpson ct al., 1995a ). 

S /H abundances are plotted in Fig. |l9L No meaningful 
Galactic gradient is seen for S/H and a very large disper- 
sion is found at any given Real, unlike the smaller dis- 
persion in the Ar/H or Ne/H Galactic abundances. The 
absolute abundances are also improbably low compared 
to the solar and interstellar abundances - on average 2 



Fig. 20. Comparison of the derived gradient slopes for 
N/O, Ne/H and Ar/H with previous determinations using 
H II regions - o ptical ( Shaver et al., 1985), IR AS ( [Simpson 
fc Rubin, 1990| ), KAO QSimpson et al., 1995b|) -, dis k plan- 



times lower that the IS sulphur abundance by Snow & 



Witt (L996); a factor of 5 for the extreme cases -, since 
sulphur is not expected to be strongly depleted into grains. 
However, the quoted sulphur abundances do not include 
the systematic uncertainties due to the effect of electron 
temperature and density (cf. Sect. |5|). In view of the rms 
densities derived from the radio observations (we note that 
the sources with the lower S /H abundances also have very 
high n c ^ ms > 10 4 cm -3 ) and the T c gradient in the Galaxy, 
we estimate that the S abundances are underestimated by 
up to a factor 4. 



7.4. Comparison of gradients with previous studies 

Fig. |2(] compares the gradient slopes derived for N/O 
Ne/H and Ar/H with previous determinations using Hn 
regions - optical (Shaver et al., 1983), IRAS (Simpson & 
Rubin, 1990| ), KAO flSimpson et al., 1995b|) -, dis k plan- 
etary nebulae of type II (Maciel fc Quireza, 1999 ) and B 
stars ( Smartt et al., 2001 ). 

Good agreement is found for N/O and Ar/H. Our 
Ne/H gradient matches the one found from planetary neb- 
ula, but is flatter than the previous infrared determina- 



etary nebulae of type II (Maciel fc Quireza, 1999| ) and B 
stars ( (Smartt ct al., 200l| ). The open stars indicate the 
slopes obtained when a Galactic T c gradient of ~ 330 
Kkpc -1 is assumed. 



tions using IRAS and KAO. Both IRAS and KAO de- 
terminations agree, however. We note that both relied on 
radio data to derive the abundances and therefore, did 
not suffer much from the systematic electron temperature 
uncertainty, while our gradients have been probably de- 
pressed because of the assumption of a single T o =7500 K. 
Indeed, if we correct our derived Ne/H gradient from the 
T e gradient in the Galaxy (the open stars in Fig. ^Oj), then 
we are in reasonable agreement with both IRAS and KAO. 

8. Summary and conclusions 

Based on the ISO spectral catalogue of compact H n re- 
gions by Peeters et al. (2001), a first analysis of the hy- 
drogen recombination and atomic fine-structure lines orig- 
inated in the ionized gas has been presented. The main 
results of this study are: 

• The SWS H I recombination lines between 2 and 8 /im 
have been used to estimate the extinction law at these 
wavelengths for 14 Hn regions. An extinction in the K 
band between and ~ 3 mag. has been derived. 

• The electron densities determined by the [O III] 
88/52 /im line ratio, in between ~ 100 and 3000 cm~ 3 , 
are shown to be substantially smaller than the rms densi- 
ties resulting from radio continuum observations. This fact 
suggests that the ultracompact cores observed in the ra- 
dio have a region of less dense gas (traced by [Om] 88/52 
/j,m) at high excitation, which could be produced by leak- 
age of radiation from the ultracompact region. This low 
density gas might be identified with the physically related, 
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Table 9. Ionic and total abundances of neon, sulphur and argon derived in the low density limit and assuming a 
single T o =7500 K. Systematic uncertainties may affect these abundances because of the T c Galactic gradient. Sulphur 
abundances for some sources could be underestimated by up to a factor 4 due to the T e and n e uncertainties (see 
Sect. |7^). 



NEON SULPHUR ARGON 

Source No+/H+ No+ + /H+ Ne/H S + + /H+ S +3 /H+ S/H Ar+/H+ Ar+ + /H+ Ar/H 

xio~ 4 xicr 4 xicr 4 xicr 6 xio~ 6 xio~ 6 xicr 6 xio~ 6 xio~ 6 



IR 01045 


0.6 ±0.3 


0.16 ±0.09 


0.7 ±0.4 


2. ± 1. 


< 0.1 


2. ± 1. 


< 0.6 


< 0.5 


< 1. 


IR 02219 


0.9±0.2 


0.9±0.2 


1.8 ± 0.3 


7. ± 1. 


0.90 ±0.09 


8. ± 1. 


0.36 ±0.03 


2.1 ± 0.2 


2.5±0.2 


IR 02575 


<0.1 


< 0.04 


<0.1 


< 0.5 


< 0.09 


< 0.6 


< 0.5 


< 0.3 


< 0.8 


IR 10589 


1.2 ± 0.2 


0.16±0.03 


1.4 ± 0.2 


9. ± 2. 


0.18 ±0.03 


10. ± 2. 


1.7 ± 0.1 


1.5 ± 0.2 


3.2±0.3 


IR 11143 


0.54±0.10 


1.1 ±0.2 


1.6 ±0.3 


12. ± 2. 


2.6 ±0.3 


14. ± 2. 


<0.2 


2.1 ± 0.2 


2.1 ±0.2 


IR 12063 


0.7±0.1 


0.59±0.10 


1.3±0.2 


5.8 ±1.0 


0.92 ±0.09 


7. ±1. 


0.43 ±0.05 


2.0 ±0.2 


2.4 ±0.3 


IR 12073 


0.34 ±0.06 


0.9±0.2 


1.2±0.2 


5.4 ± 1.0 


1.9 ±0.2 


7. ±1. 


0.12 ±0.04 


1.9 ±0.2 


2.0±0.2 


IR 12331 


1.9±0.3 


0.59±0.10 


2.4 ±0.4 


19. ±3. 


0.53 ±0.05 


19. ±3. 


1.6 ±0.2 


3.3 ±0.4 


4.9 ±0.6 


IR 15384 


1.8±0.3 


0.25 ±0.05 


2.0 ±0.4 


9. ±2. 


0.19±0.02 


9. ±2. 


2.6 ±0.2 


1.8 ±0.3 


4.4 ±0.5 


IR 15502 


1.2±0.3 


0.13±0.03 


1.4 ±0.3 


3.4±0.7 


0.033 ±0.006 


3.4±0.7 


2.7 ± 0.4 


0.6 ±0.1 


3.4±0.5 


IR 17221 


2.4±0.5 


0.04 ±0.01 


2.4 ±0.5 


12. ± 2. 


0.044 ±0.010 


12. ± 2. 


4.4 ±0.6 


1.2 ± 0.2 


5.6±0.8 


SgrC 


3.0±0.7 


<0.05 


3.0±0.7 


7. ±2. 


<0.08 


7. ±2. 


11. ± 2. 


<0.4 


11. ± 2. 


IR 17455 


1.6±0.3 


0.31 ±0.05 


1.9 ±0.3 


12. ± 2. 


0.19 ±0.03 


12. ± 2. 


2.7 ±0.2 


1.8 ±0.2 


4.4 ±0.4 


IR 18032 


2.0±0.4 


<0.02 


2.0±0.4 


8. ±1. 


<0.03 


8. ±1. 


5.3 ±0.6 


0.71 ±0.09 


6.0±0.7 


IR 18116 


2.3±0.4 


0.08 ±0.02 


2.4 ±0.4 


15. ±3. 


<0.03 


15. ±3. 


7.4 ±0.6 


0.9 ±0.1 


8.2±0.7 


IR 18317 


2.6±0.5 


0.08 ±0.01 


2.6±0.5 


13. ±2. 


0.034 ±0.009 


13. ±2. 


4.9 ±0.4 


1.1 ±0.1 


6.0±0.5 


IR 18434 


2.2±0.4 


0.26 ±0.05 


2.4 ±0.4 


7. ±1. 


0.13±0.02 


7. ±1. 


3.0 ±0.3 


1.7 ±0.2 


4.6±0.5 


IR 18479 


1.5±0.3 


0.30 ±0.06 


1.8 ±0.3 


3.7±0.7 


<0.05 


3.7±0.7 


1.7 ±0.2 


<0.2 


1.7±0.2 


IR 18502 


2.6±0.4 


0.26 ±0.05 


2.8±0.5 


10. ±2. 


0.12 ±0.02 


10. ±2. 


2.9 ±0.3 


1.1 ±0.1 


3.9±0.4 


IR 19442 


2.0±0.9 


<0.03 


2.0±0.9 


4. ± 2. 


<0.1 


4. ± 2. 


6. ±2. 


<0.3 


6. ±2. 


IR 19598 


0.37 ±0.06 


0.38 ±0.06 


0.8 ±0.1 


1.7±0.3 


0.15 ±0.01 


1.8±0.3 


0.43 ± 0.03 


0.61 ±0.06 


1.05±0.10 


DR21 


0.9±0.2 


0.11 ±0.05 


1.0±0.2 


2.6±0.5 


<0.02 


2.6±0.5 


2.9 ±0.4 


0.35 ±0.07 


3.2±0.5 


IR 21190 


0.7±0.1 


0.29 ±0.05 


1.0±0.2 


3.0 ±0.5 


0.20 ±0.02 


3.2±0.5 


0.50 ±0.06 


1.4 ±0.1 


1.9±0.2 


IR 22308 


0.9±0.2 


<0.02 


0.9±0.2 


5.5 ±1.0 


<0.05 


5.5 ±1.0 


2.8 ±0.5 


0.6 ±0.2 


3.5±0.7 


IR 23030 


1.4±0.2 


0.09 ±0.01 


1.5±0.2 


10. ±2. 


0.17±0.02 


11. ± 2. 


1.0 ±0.1 


2.0 ±0.2 


3.0±0.3 


IR 23133 


1.1 ±0.2 


<0.01 


1.1 ±0.2 


6. ±1. 


0.025 ±0.006 


6. ±1. 


2.9 ±0.3 


0.73 ±0.08 


3.7±0.3 



low brightness envelopes recently mapped in radio around 



ultracompact Hn regions (Kurtz et al., 199£; Kim & Koo. 
200lj)~| 

• Significant variations are found in the degree of ion- 
ization of the H II regions as measured by the ratios 
Ar++/Ar+, N++/N+, S+ 3 /S++ and Ne++/Ne+. These 
ratios, which span a range of ionization potential up to 
41 eV, correlate well with each other, suggesting that the 
spectr al hardening affects equally the full range of ionizing 



of radiation pressure effects during the collapse phase (cf. 
Kahn, 1974; Shields fc Tinsley, 1976| ). Last, the hardening 
of the radiation field with Real may reflect a systematic 
variation in the "age" of the star powering the H n regions. 
Possibly, in the inner parts of the Galaxy ionizing stars 
evolve more "rapidly" (ie. are cooler for the same luminos- 
ity) than in the outer Galaxy. In that respect, the detailed 
photoionized model for the H II region IRAS 18434-0242 
(G29.96-0.02), using recent line-blanketed non-LTE stel- 
lar atmosph ere models including winds ( Schaerer fc de 



energi es. In addition, the ionic ratios vary with Rn,i. with 
the Hn regions in the outer Galaxy being powered by a 
seemingly harder radiation field. This increase in degree 
of ionization, ie. in properties of the ionizing star, could 
be related to the metallicity gradient. However, the pre- 
cise mechanism coupling the stellar characteristics to the 
elemental abundances remains unclear. Elemental abun- 
dances may couple directly to the stellar atmosphere opac- 
ity which controls the flux and its dependence on the 
frequency. Qualitatively, theoretical models support this 
coupling (e.g. Balick & Sncdcn, 1976). The opacity by 



Koter, 1997 ), finds good agreement with all the ion pair 



trace ions also controls the stellar wind, which in their 
turn have a major influence on the emerging stellar flux. 
Again, models sho w qualitative agreement (e.g. |Schaerer 
& de Koter, 1997j ). Third, the connection between ele- 
mental abundances and stellar flux may be indirect. For 
example, the elemental abundance may control the dust 
opacity and thereby the mass of the star formed because 



ratios using a luminous star which has already evolved off 
the main sequence (Paper III). In such a model, the varia- 
tions in the degree of ionization with Real might reflect a 
selection effect in our sample where, for example, the rate 
of expansion of the compact Hn region is linked to the 
environment (ie. metallicity). 

• A gradient in N ++ /0 ++ is observed in the sense of 
a decreasing ionic abundance ratio with Real- This ob- 
served ionic abundance ratio N ++ /0 ++ has been con- 
verted into the elemental abundance ratio N/O using the 
observed ionization ratio Ne ++ /Ne + , sensitive to the ion- 
ization state of the H n region, and detailed photoioniza- 
tion model calculations. The N/O abundance ratio also 
shows a decrease with Real, although shallower than in 
the case of N ++ /0 ++ , pointing to the existence of an in- 
trinsic variation in the gas phase N/O abundance. Thus, 
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there seems to be an important contribution to the syn- 
thesis of nitrogen by secondary production. The derived 
slope for N/O (AlogN/O = -0.056 ± 0.009 de xkpc- 1 ) 
is in agreement with previous H H regions (e.g. [Bimpson 



et al., )_995b ) and B stars (e.g. smartt et al., 2001) deter- 
minations. 

• Our analysis method underestimated the sulphur abun- 
dances by up to a factor 4 and hence, no meaningful 
Galactic gradient could be derived. We estimated a lower 
limit for the oxygen abundance through the derivation of 
++ /H + . We miss the important On ionization stage, 
which does not have any infrared fine-structure line. On 
the other hand, because the [Nil] 122 /im line is very 
weak and it was only detected in a limited number of 
sources, and no radio information was available for all 
these sources, N/H could only be determined for a few 
sources and therefore, a gradient could not be drawn. We 
found, however, supersolar nitrogen abundances for the 
sources in the inner 5 kpc. 

• The total elemental abundances of neon and argon 
with respect to H have been also derived. Ne/H and 
Ar/H have been found to decrease with the galactocen- 
tric distance as AlogNe/H = -0.039 ± 0.007 dexkpc" 1 
and AlogAr/H = -0.045 ± 0.011 dexkpc" 1 . The slopes of 
neon and argon coincide within the uncertainties, indicat- 
ing that both species have a similar stellar nucleosynthe- 
sis history. However, the existence of a Galactocentric T e 
gradient may enlarge these elemental gradients. Because 
only a few radio recombination line observations towards 
the Hn regions, which can be used to characterize the 
nebular T c , are available, we do not have enough infor- 
mation to correct individually the derived abundances. 
However, adopting a decrease in T c of ~ 330 K kpc -1 
with galactocentric distance, we estimate that the slopes 
for the Ne/H and Ar/H gradients become —0.06 and —0.07 
dexkpc -1 , respectively. These corrected slopes are in rea- 
sonable agreement with the previous determinations by 



IRAS ( Simpson & Rubin, 199C| ) and KAO ( |Simpson et al 



1995b[j. Ar/Ne, which is independent of T e , is found to be 
approximately solar and constant with Real- 

The complete infrared spectra from 2.3 to 196 /im pro- 



vided by the ISO catalogue of galactic H n regions (Peeters 
et al., i 001) has allowed us to study the ionizing conditions 



in these sources and derive relative and absolute elemen- 
tal abundances across the Galactic disk. The first analysis 
presented in this paper has shown that useful conclusions 
can be drawn from these data. However, there are still 
open issues and to improve on the present results addi- 
tional information is required. In particular, the detailed 
knowledge of the extent and relative spatial distributions 
of the fine-structure atomic lines is lacking for almost all 
the compact Hn regions. The knowledge of the electron 
temperature is also critical in refining the elemental abun- 
dances. Observations of radio recombinations lines will be 
helpful to estimate T D for all the sources discussed in this 
paper and to correct for the T e dependence. 
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